Certifying entangled measurements in known Hilbert spaces 
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We study under which conditions it is possible to assert that a joint demolition measurement can- 
not be simulated by Local Operations and Classical Communication. More concretely, we consider 
a scenario where two parties, Alice and Bob, send each an unknown state to a third party, Charlie, 
who in turn interacts with the states in some undisclosed way and then announces an outcome. We 
show that, under the assumption that Alice and Bob know the dimensionality of their systems, there 
exist situations where the statistics of the outcomes reveals the nature of Charlie's measurement. 
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I. INTRODUCTION 

Quantum nonlocality, the fact that entangled quantum 
systems separated in space can violate Bell inequalities 
is one of the most profound discoveries in science. 
This nonlocal nature of entanglement has also been iden- 
tified as an essential resource for various quantum infor- 
mation processing tasks. It allows to reduce communica- 
tion complexity or makes it possible to devise device- 
independent protocols for quantum key distribution [||, 
genuine random number generations Q , and state tomog- 
raphy Q. Indeed, these tasks can be performed without 
resorting to the actual inner working of the devices, which 
is a feat without counterpart in the classical world. 

In all of the above instances, quantum nonlocality 
arises when space-like separated measurements are per- 
formed over jointly prepared entangled states. However, 
the existence of nonlocality can also be manifested in a 
dual setup [^,[0: quantum states which arc prepared 
separately (hence unentangled) exhibit anomalous cor- 
relations when measured jointly. In this case, more in- 
formation can be revealed by joint measurements than 
can be gained by using any sequence of Local Operations 
assisted by Classical Communication (LOCC). In partic- 
ular, the superiority of globally entangled measurements 
over LOCC ones has been proved conclusively The 
improved performance is due to entanglement occurring 
in the process of measurement and again it is not associ- 
ated with the states. 

For quantifying the effectiveness of entangled measure- 
ments, different figures of merit have been proposed, like 
the Shannon mutual information |^ or the quantum fi- 
delity [||. A common feature of those is, though, that one 
must rely both on the precise form of the states prepared 
and on the a priori probabilities of the preparations; oth- 
erwise, the estimations on the figure of merit would be 
not reliable. Hence, in order to ascertain that a measure- 
ment is truly globally entangled and not only an LOCC 
one, one must resort to a detailed knowledge of the prepa- 
ration procedure. Clearly, this issue sheds doubt on the 
reliability of the obtained results. 

In the present paper we take another approach, more 
in the spirit of a black box scenario, and tackle the prob- 
lem of certifying entangled measurements by introducing 
a kind of dual Bell inequalities. These inequalities involve 



correlation terms which can be gathered from experimen- 
tal data, and do not depend on the specific form of the 
quantum states to be prepared. Actually, the only nec- 
essary condition entering in the derivation is the Hilbert 
space dimension of the prepared particles. Taking into 
account that entangled measurements enable increased 
classical capacity of quantum channels |^ and efficient 
quantum state estimation JTo| , in the future it would also 
be interesting to find applications related to these tasks 
within a black box approach. 

The notation we use along the article is introduced 
in the next section. Section [II provides an instance of 
certifiable entangled measurements based on a nonlinear 
inequality. Section |^ is devoted to a numerical study 
where linear inequalities are applied to detect entangled 
measurements in the simplest possible scenarios. Sec- 
tion 1^ summarizes our conclusions. 



II. PRELIMINARIES AND NOTATION 

Let us imagine the following communication scenario, 
close to the spirit of the simultaneous message passing 
model [0. Two separated parties, Alice and Bob, have 
each some preparation device with N possible inputs. 
Alice (Bob) is thus able to prepare any of the unknown 



qudit states {p^ G C"}x=i iWy G <C'^}y=i)- Alice and 
Bob then send their states to a third party, Charlie, who 
performs a measurement over those two states, announc- 
ing a dit a of dimension K at the end of the process. 
Denote by {Ma : a = 0, K — 1} the elements of Char- 
lie's Positive Operator Valued Measure (POVM). Then, 
for each pair of inputs x and y, several repetitions of 
this primitive would allow Alice and Bob to estimate the 
frequency of occurrence 

P{a\x,y) = tr{p^(E)(TyMa) (1) 

of each outcome a = 1, 2, . . . , A". Depending on the form 
of the POVM elements {Ma : a = 0, K - 1}, we will 
distinguish four different types of measurements: 

Classical measurements. Suppose that Charlie's mea- 
surement is fixed in advance to a given complete pro- 
jective measurement on each system (say, in the com- 
putational basis), followed by some random processing 
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FIG. 1: A schematic picture of the 2-dimensional sUce of 
the probabihty space for different sets in the general NDK 
scenario. The square represents the classical set , the circle 
and the ellipse designate the set of unentangled (p""*!"') and 
the set of entangled measurements (P''"*), respectively. In 
the purple region entangled measurements are witnessed by 
linear inequalities. 



of the raw output. Any set of distributions P{a\x,y) 
so generated can be reproduced in a classical setting 
where Alice and Bob send each one dit of information 
to Charlie, who, in turn, applies a random function 
/ : {0, ...,D — 1} X {0, D— 1} — >■ K and then announces 
the result. For each NDK scenario, we will denote by 
the set of all possible correlations P{a\x, y) that can 
be attained in such a way. 

It is easy to see that any point of P*^ is a convex combi- 
nation of deterministic classical strategies. As we will see, 
though, the converse is not true. On the other hand, if 
Alice, Bob and Charlie share random variables, then any 
convex sum of deterministic strategies can be attained. 
Since for any NDK scenario the number of deterministic 
classical strategies is finite, the set of correlations attain- 
able through classical maps and shared randomness is a 
polytope, that we will denote by P*^. This classical set 
is symbolized in Figure |^ by a square. 

Unentangled and LOCC measurements. If each of the 
POVM elements Ma is a separable operator, wc will 
say that the measurement M is unentangled. We will 
say that Charlie's measurement can be attained via Lo- 
cal Operations and Classical Communications (LOCC) 
if M corresponds to a sequence of local measurements 
on Alice's and Bob's individual qubits, with each mea- 
surement possibly depending on the outcomes of earlier 
measurements. It is known that the class of unentangled 
measurements is strictly greater than the class of LOCC 
measurements Q. This implies that P""'^"* D pLOCC ^ 



where p^^^^^t ^ pLOCC j-^^p denote the sets of distri- 
butions P{a\x, y) attainable through unentangled and 
LOCC measurements in a given NDK scenario. It is 
an open question whether such an inclusion is strict. 

As before, if we allow Alice, Bob and Charlie share 
some prior random variables, the resulting sets of cor- 
relations P""<="* D pLOCC g^j-g convexifications of the 
former. A two-dimensional slice of the unentangled set 
punent jg symbolizcd in Figure |^ by a circle. 

General measurements. Here the measurement oper- 
ators are only limited by positivity and normalization 
constraints, i.e., 

K-l 

Ma to, J2Ma^l, (2) 

a=0 

where the operators Ma G B{C^ (g) C^) may very well be 
entangled. These are the most general measurements al- 
lowed by quantum mechanics. In analogy with the former 
sets, we will denote by P*^"* and P'^"* the set of distribu- 
tions accessible through general joint measurements and 
its convex hull. The latter is depicted as an ellipse in 
Figure |l[ 

The numerical study of the sets P*^ , P""ent ^^^^ pent -g 
much more convoluted than that of their convexifications 
pC ^ punent ^^^^ pent ^ Consequently, most of this article 
is devoted to the numerical analysis of the latter. 

If we reflex a bit about the previous definitions, it will 
soon be clear that the above sets of correlations can only 
be distinguished experimentally in a black box way if 
the NDK scenario is such that N > D. Otherwise, by 
forcing Alice and Bob to send orthogonal states, Char- 
lie could always infer the values oi x,y through an ap- 
propriate projective measurement, and so any conceiv- 
able distribution P{a\x,y) could be classically realized. 
The simplest scenario where we may expect to find an 
interesting structure is thus the 322 scenario. And, in- 
deed, the next example shows that there the sets P*^ and 
punent (jjffgj. froni p<="*. ft also proves that neither of 
the former two sets is convex, and so the convexifications 
pC pLOCC punent g^^^ nou-trivial. 



III. AN INSTANCE OF ENTANGLEMENT 
DETECTION 

Consider the 322 scenario, and denote by P the matrix 
of probabilities P^y = P{0\x, y). Suppose now that 

/ 1/4 l/4\ 
P= 1/4 1/4 . (3) 
\l/4 1/4 j 

This set of probabilities can be attained if the states 
{px, <Jy} and the POVM element M = Ma=o correspond 



3 



to 



pi = |0)(0|, P2 = 

ai = |1>(1|, CJ2 = 



P3 = I I, 
<^3 = I I, 



(4) 



(5) 



where = ^TjdOO) + |11)). It hence foUows that 

p g pent 

We will prove that Af must be an entangling mea- 
sure (i.e., P ^ punent-^ reducUo ad ahsurdum. 
Imagine, therefore, that M is separable, i.e., AI = 
X^iLi ^i\ui){ui\®\vi){vi\. Then, the condition Pn = im- 
plies that: 1) pi or ai (or both) is a pure quantum state; 
and 2) there exists a subset of indices X C {1, 2, A'} 
such that 

A/ = ^ X,pi ® \v,){v,\ + A,|u,)(ii,| ® ai- (6) 

Here denotes I2 — o;^, for any qubit state uj. From 
tr(Mpi ® (Ti) = Pii ^ P12 = tr(Af/9i (g) 02) we have 
that o"! 7^ (72, and thus tr(cri(T^) 7^ 0. Likewise, 
tr(px/oi>),tr((Tj,(T,^) 7^ 0, for aU y ^ y',x ^ x' . This 
implies that 



tr(A//p2 ® ^12) = ^ A,tr(pjLp2)tr(|v,)(w,|cr2) 
+ ^ Ai tr ( I )(u J I /CI2 ) tr ( (T J*- (72 ) 



> (5^A,tr(|u,)(w,|a2) +(5^A,tr(|w,)(w,|p2), (7) 

for some 5 > 0. The condition P22 = therefore requires 
that tr(|wi)(vj|a2) = 0,i e I, and tr(|Mi)(ui|p2) = 0,i ^ I, 
so 

A/ = ^A,/9]^«)o-^ + ^Aj/9^(8)cr]^. (8) 

Using this last decomposition, we get 

tr(Afp3 ® (73) =^Aitr(p]^p3)tr(cr^ era) 
■I 

X! '^»tl'(P2"/'3)tr(cr]^cr3) 



i6l 



if 



> (J'Xl'*^* - ® ^^2) > 0. (9) 

with 5' > 0. But P33 = 0, from which it follows that M 
is indeed entangling. Actually, it is possible to obtain a 
non- linear witness to detect Af's entanglement if we just 
follow the previous steps carefully while keeping track of 
the approximation errors (see the Appendix) . 

We have thus proven that P ^ punent ^ j ^j^g matrix 
of probabilities P cannot be observed if only unentangled 



measurements are allowed. However, the matrix P^y can 
be expressed as a convex combination of four determin- 
istic classical strategies. Indeed, 



P 



1 



^0 1 1 

.0 




^0 0> 




^0 0; 
(10) 

The sets pc ^ pLOCC ^ punent g^j.^ ^^^^^ jjQ^ convex, 
since otherwise P would belong to them. Allowing the 
three parties to have shared randomness thus leads to a 
radically different scenario. 



IV. 



NUMERICAL STUDY OF , p""™' AND 
pent jjvj rpjjg SCENARIO 



A. Procedure 

Since convex sets can be completely characterized by 
systems of linear inequalities, it is legitimate to investi- 
gate how different linear constraints on P'^ of the form 
y ^xyPxy < Wc niay be violated by elements of 

punent gj^f^ pent 

This leads us to the problem of maximizing quantities 
Uke 



N 

E 

x,y=l 



^^x ,y Px ^y 



N 

E 

x,y=l 



W^ytr{p^(g,ayM) (11) 



over all possible POVM elements M in the unentangled 
or the general class, and over all qubit states p^ and 
ay so as to get the values Wunent and Went- Clearly, if 
Went > Wunent for some matrix W, then there exist exper- 
imental situations where one can prove that the statistics 
observed in the lab cannot be simulated with unentan- 
gled measurements and shared randomness. 

Numerical optimization to obtain Wunent is carried out 
very similarly to the iterative algorithm used in ||l^ . The 
steps are the following: 

1. Generate some random pure qubit states px, (Jy, 
x,y=l,2,...,N. 

2. In dimension Z) = 2, for fixed states px,cry, max- 
imizing Eq. ( pi] ) reduces to a semidefinite prob- 
lem: define thus F = "Y^^ y^i^x.yPx ® cFy, and 
maximize tr(AfF), subject to Af ^ 0, 1 - A/ ^ 
0, PT{M) ^ 0, 1 - Pt(m\^ 0, where PT denotes 
Partial Transposition 
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3. Fix M and ay and maximize tT{GxPx) = 
{ipx\Gx\ipx) for X = l,2,...,iV where Gx = 
trs(X]y 1 ^ ayM). This amounts to find the max- 
imum eigenvalue of Gx with the corresponding 
eigenvector \i}}x) for each x. 

4. For fixed M and p^, the optimal ay can be found 
such as px in step 3. 
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5. If convergence has not been achieved, go back to 
step 2. 

This optimization algorithm may encounter several lo- 
cal optima, therefore we must iterate it many times in 
order to ascertain with reasonable confidence that the 
largest maximum has been found. In Sections [VB| and 
IV C the semidefinite programs described in step 2 were 
solved using the SeDuMi package . 

A very similar procedure can be applied if we wish 
to inquire about the maximal value Went which can 
be achieved within the class of general quantum mea- 
surements. The only difference is in step 2, where 
we have to maximize tr(Mi^) subject to the only con- 
straint ^ il/ ^ 1. However, in this case the max- 
imum for a fixed F can be easily found by diagonal- 
ization F ~ resulting in the optimal 

M = Eti ^^^^^\(l>i){cl),\- Thereby the optimization 



algorithm for the case of entangled measurements is much 
faster. 



B. The 322 scenario 

We next describe the way linear witnesses were pro- 
duced for entangled measurements. Our starting point is 
the classical polytope P'-^ for the 322 scenario, consisting 
of 104 extremal points in dimension 9. All its facets can 
be enumerated by using the double description method 
implemented in Fukuda's cdd package [p^ . The polytope 
has 1230 facets; among them, 9 describe positivity 
facets Px.y > 0, x,y = 1,2,3. The rest of them define 
non-trivial inequalities W. However, most of them are 
equivalent under the following operations, 

• permuting Alice's inputs x 

• permuting Bob's inputs y 

• exchanging parties Alice and Bob 

• multiplying by —1 all coefficients Wx,y of the ma- 
trix W. 

It turns out that there are 13 inequivalent facets of 
the polytope P*-^. Those are listed in Table |, where 
we also give the maximum classical value Wc, the value 
Wunent achievable with unentangled measurements, the 
value Went attainable with general measurements, and 
the number of vertices lying on the particular facet. We 
must stress that the values for Wunent and Went come 
from numerical optimization, and so only local optimal- 
ity is guaranteed. Nevertheless, due to the small dimen- 
sionality of the problem and the number of iterations of 
the main algorithm, we are quite confident about their 
overall optimality as well. Interestingly, we found that in 
each of the 13 cases Wunent could be achieved with rank-2 
projective measurement operators. 

We now focus on inequality #4 from Table ||, which is 
the only inequality for which Went > wunent, and hence 



it enables to witness entangled measurements. The in- 
equality looks (in an equivalent form) as 



— Pll — Pl2 

^ 2-f 3^/6 ^ 



- Pl3 + P2I 

2.3371, 



P 



23 



Pii — P?.2 — P: 



33 



(12) 



and can be violated by entangled measurements up to the 
value of 2.5. Next, we present the actual measurements 
and states achieving Went and Wunent- 

Entangled case. The POVM clement M is a rank-2 
projector 



M = \m){m\ + \m-^)(m^ 



with 



(0, cos 61, sin 61,0) 
cosd . , . 



V2 



, — sm (p sm W, sm (p cos f 



cos <p 



(13) 



(14) 



having — 2 arctan( ^ ) and — n/3. The partial 
transpose of M has the following eigenvalues: 



A 



5- \/4T 14 5 
10 ' 5' 5'~ 



10 



(15) 



with the first entry being negative. Both Alice's and 
Bob's states are pure and real valued, px = \'4'x){4'x \ and 
ay = \4'v){4'y\i respectively. Their explicit forms are, 
respectively. 



1^1) 
1^2) 

IV^3) 



and \(l)y) 



= (1,0), 

: (cos a, sin a), 
: (cos 2a, sin 2a), 

for y = 1,2,3 with a = 
± 0.9117 rad. With these val- 
ues in hand, the probability matrix P with components 
Px,y = tr(/9:r ^ PyM) looks as follows: 

^0.125 0.25 0.875\ 

0.25 0.5 0.75 . (16) 
*^0.875 0.25 0.125/ 



= \^y 

2arctan(2vl0p^ 



P = 



And, indeed, evaluating the left-hand side of inequal- 
ity (y|), we get Went =2. 5. 

Unentangled case. In appropriate local bases, a sep- 
arable rank-2 projective operator can be written in the 
form 



where |m 
i arcc 
and (7. 



M = |0)(0| ® \m){m\ + |1)(1| ® |0)(0|, (17) 

) = (cos 6, sin 6), and in our particular case, 9 = 
~ 0.911738 rad. The states = \ipx){'ipx\ 
\(l)y){(py\, on the other hand, are given by 

(1,0), 



l^l> 


= (1,0), 




H2) 


= (1,0), 


\h 


1^-3) 


= (0,1), 





(cos a, sin a), 
(cos —a, sin —a). 



(18) 
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TABLE I: Results for the 13 facet inducing inequalities in the 322-scenario. 



Case 




'^unent 


Went 


vertices 




Wi2 


Wi3 


W21 


W22 


^^23 


W31 


H/32 


W33 


1 


6 


6.4006 


6.4006 


11 


-2 


-1 


1 


2 


4 


-2 


4 


-5 


-1 


2 


2 


2.8284 


2.8284 


13 


-2 





2 


-2 


1 


-1 





1 


-1 


3 


1 


1.4142 


1.4142 


20 


-1 


-1 








1 





1 


-1 





4 


2 


2.3371 


2.5 


16 


-1 


-1 





-1 


1 


1 


1 


-1 


1 


5 


1 


1.3371 


1.3371 


12 


-1 


-1 


1 


-1 





-1 





1 





6 


2 


2.6742 


2.6742 


10 


-2 


-1 


2 


-2 


1 


-1 





1 





7 


2 


2.1623 


2.1623 


9 


-1 





1 





2 


-2 


1 


-2 


-1 


8 


2 


2.1403 


2.1403 


10 


-2 


-1 


1 





2 


-2 


2 


-3 


-1 


9 


1 


1.2058 


1.2058 


10 


-1 


-1 


1 





1 


-1 


1 


-2 


-1 


10 


2 


2.7275 


2.7275 


14 


-3 


-2 


1 


-2 


2 





1 


-1 


1 


11 


1 


1.3094 


1.3094 


16 


-2 


-2 


1 


-1 


1 





1 


-1 





12 


3 


3.8467 


3.8467 


11 


-5 


-3 


2 


-3 


3 


-1 


2 


-1 


1 


13 


2 


2.1186 


2.1186 


9 


-3 


-1 


2 


-1 


2 


-3 


2 


-3 


-1 



with 



- arctan I 



-1.11493 rad. With 



these, the probabihty vector is 



'0.375 P12 Pi3 
P = I 0.375 P12 Pi3 
1 P12 P12 1 



(19) 



where Pi 2 



0.193814 and P13 = ^^(5 



4(4+%/6) 

2-\/6) — 0.959279. Plugging these numbers into the left- 
hand side of inequality (|l2[), we get Wunent = — 
2.3371. In some sense, the probability vector ( [l9|) ap- 
proximates the best way among unentangled measure- 
ments the probability vector (Uw corresponding to an 
entangled measurement. 



C. The 422 scenario 

The classical polytopc P^ for the 422 scenario is 
spanned by 520 non-redundant vertices in the 16 dimen- 
sional probability space. However, it turned out that the 
full characterization of this polytope in terms of facets 
is computationally an elusive task. Instead, the prob- 
lem was approached in a different way. We scanned 
through all the inequalities with small integer coefficients 
(—1,0,-1-1), and sorted out all of them which are facet 
defining. Several inequivalent facets have been found in 
this way (beyond the ones, which arc just straightfor- 
ward extensions of the 322-type inequalities in Table |). 
In Table |l[ we list all those 10 inequivalent facets for 
which Went > Wunent, licnce witncssiug entangled mea- 
surements. Among them, #9 is equivalent to the single 
322-type witness in (p^). Table || contains results on 
Wc, Wunent, Went, and also glvcs the sole negative eigen- 
value Ai of PT{M) achieving Went- Since all the other 
eigenvalues are positive, the absolute value |Ai| defines 
negativity, a valid entanglement measure U%- 



Inequality #2 is interesting on its own. In this case. 
Went = 3.3195 is attained with a negativity of 0.1744, 
although a larger value of negativity of 0.1777 is found 
at another locally optimal point with w^^^ = 3.145. This 
feature resembles the dual case obtained in a Bell sce- 
nario, where for several tight 2-party Bell inequalities 
the maximum quantum value was achieved with non- 
maximally entangled 2-qubit states [|l8| . It has also been 
observed that for all the 422 witnesses for which unentan- 
gled measurements gave numerically maximal violation, 
the measurements could always be brought to a form of 
rank- 2 projective matrices. 

To conclude this section, the 422 scenario gave us a 
few new entangled measurement witnesses over the sin- 
gle witness of ( p^ of the 322 scenario. However, no sig- 
nificant improvement could be found in the efficiency of 
detecting entangled measurements with respect to unen- 
tangled measurements. Also, no simple generalization of 
inequality ( p2[ ) for larger input alphabets was discovered 
among the 422 inequalities. 



V. CONCLUSION 

In this paper, we have studied how to translate the con- 
cept of witnessing entangled measurements to the device- 
independent arena. We found that even in a black-box 
scenario it is possible for a theorist to assert that a joint 
demolition measurement cannot be simulated with any 
sequence of local measurements assisted by classical com- 
munication: the only assumption one has to rely on is 
the dimensionality of the probe states. Following this 
line of thought, we derived correlation inequalities which 
allow to test the superiority of entangled measurements 
over u nentangl ed ones in two-qubit experiment s whe re 9 
(Sees. HI IV B) and 16 correlation terms (Sec. IVC) are 
estimated. 

Previous routes to certify the 'non-locality' of unknown 
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TABLE II: Results for the 10 facet inducing inequalities witnessing entangled measurements in the 422-scenario. 



Case 




'^unent 


Went 


vertices 


Ai 




Wi2 


Wi3 


Wi4 


W21 


W22 


W23 


W2i 


W31 


W32 


W33 


W34 


W41 




W43 


W^44 


1 


2 


2.3371 


2.3510 


31 


-0.1238 


-1 


-1 


-1 





-1 





1 





-1 


1 


-1 


1 


1 


-1 





1 


2 


3 


3.2361 


3.3195 


34 


-0.1744 


-1 


-1 


-1 





-1 





1 


1 





1 


-1 


1 


1 


-1 





1 


3 


2 


2.3371 


2.4369 


38 


-0.1154 


-1 


-1 


-1 


1 


-1 





1 


-1 

















1 





1 


4 


2 


2.3371 


2.4339 


25 


-0.1133 


-1 


-1 


-1 


1 


-1 





1 


-1 








1 


1 


1 


-1 





-1 


5 


1 


1.3371 


1.3413 


28 


-0.1016 


-1 


-1 








-1 





-1 














1 





1 


-1 


-1 


6 


2 


2.3371 


2.4322 


24 


-0.1089 


-1 


-1 








-1 





-1 


1 





1 


-1 


-1 


1 








1 


7 


2 


2.1623 


2.3028 


19 


-0.0718 


-1 


-1 








-1 


1 


-1 


1 


1 





-1 





1 








-1 


8 


2 


2.1623 


2.3028 


26 


-0.0718 


-1 


-1 








-1 


1 








1 


-1 


-1 


1 


1 


1 


-1 


-1 


9 


2 


2.3371 


2.5 


72 


-0.1403 


-1 


-1 








-1 


1 





1 














1 


-1 





1 


10 


2 


2.2058 


2.3773 


22 


-0.0982 


-1 


-1 








-1 


1 





1 





1 


-1 


-1 


1 


-1 


-1 


1 



measuring devices involved carrying out complete tomog- 
raphy of their associated POVM elements or en- 
visaging a state discrimination problem where entangled 
measurements have an advantage with respect to unen- 
tangled ones Either scheme has been conducted 
experimentally ^ with success. Both approaches, 
though, relied crucially on a detailed knowledge of the 
probe states to be measured, and so their conclusions 
cannot be considered definite. 

In contrast, in this article we have proven that the ex- 
perimental inaccuracies arising from an imperfect prepa- 
ration of the probe states can be completely eliminated 
by certifying entangled measurements within a black box 
approach. Given the simplicity of our inequalities, we 
thus propose their experimental implementation as an 
interesting challenge. 
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Appendix A: A Non-linear witness 

In order to obtain a witness to certify the non-locality 
of P in Eq. (^, we will make use of the following result. 

Proposition 1. Let A e B{C'^) satisfy ^ A^ld, and 
let uJi,UJ2 €: B{C^) be two normalized quantum states, 
with tr{Au}i) = Pi, for i — 1,2. Then, 

tr{uji-U2)<{f{Pi,P2))\ (Al) 

where 

/(Pi,P2) = ^Pi ■ P2 + v/(l - -Pi) • (1 - P2). (A2) 

Moreover, there exist two states uji,uj2 G B{C^) and a 
POVM element A E B{C^) that saturate the former in- 
equality. 

Proof. Let fii\i){i\ be the spectral decomposition of A, 
where < fii < 1 and \\i)}fz^ is an orthonormal basis 
of C^. Then, Eq. ( |Al|) follows from the next chain of 
inequalities: 



(tr{c.i . 0:2} f = (E.,,(*l^ib)(.7V2N))'^' < (E.,, i{^M){J\c^l\J))'^' i{^\u2\^){3\u^2\J))'^'] - 

= V7W¥ + - Pi) ■ (1 - P^)- 
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(A3) 



To see that the bound is optimal, choose A ~ |0)(0| □ 
and notice that the relations |(i/'i2|0)P ~ Pi 2 and 

|(V'i|V'2)P = f{Pi,P2f hold for the states |V^1,2) = Taking A = M and cji = ® ai [tu^ = pi ® a^), 

V^^|0) + Jl - Pi 2II). ' = Pk®cri (w2 ^ pi® Uk) in Proposition ||, it is easy 

to see that 
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tT{pj ■ pk) < mmfiPjuPkif 
(tr(a,-(7fc) < minf{Pij,Pikf). 



(A4) 



These relations will play an important role when com- 
bined with the next proposition. 

Proposition 2. Let uj = \\u){u\ ® \v){v\ G ^(C^ (g) C^), 
with A > 0, let {pj, crj}j=i be normalized qubit states such 
that tr{pj ■ pk}, tr{aj ■ ak} < C, and let tr{ujpj ®<Jj} < c. 
Then, 



X < 



4:C 



(1 - ■ 

2 



(A5) 



Proof. Let pj = ^. _ i+n^o;^ hypothesis we 

have that 



{I + u ■ rfij ) ■ {1 + V ■ fij) < c, 



(A6) 



for i = 1, 2,3, where c := 4c/A. It follows that, for each 
inequality j, at least one of the factors on the left-hand 
side is smaller or equal than This implies that there 
exist two indices j, k e {1, 2, 3}, j ^ fc, such that either 



(1 4- It • TOj.fe) < v5 



(1 + t/- fij^k) < 



(A7) 



(A8) 



holds. Let us assume that Eq. (|A7|) is true. Then we 
have that 

II mj + rfikW > —u- {fhj + rhk) > 2 — 2\/5, (A9) 
and hence 

4v^ 



A < 



On the other hand. 



2 - ||too -I- mk\ 



(AlO) 



\\rhj+7nkf < 2{l + rnj-rnk) = 4tr(pj-pfc) < AC. (All) 



Combining Eqs. ( AlO ) and ( All ), we arrive at 

4c 

A < 



(1-VC) 



(A12) 



Had we supposed that Eq. ( |A8| ) wer e tr ue, we would 
have obtained the same relation. Eq. (A5) then follows 
from the fact that at least one of the conditions (|A7|), 
(|A^ ) must hold. □ 

We are now ready to derive the non-linear inequality. 
Suppose that M is of the form 
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M = Y, \\ur){u,\ ® \v,){v,\, (A13) 

i=l 

with A; > 0. It is straightforward that the relation 

3 



Xi{ui\pj\ui)(vi\aj\vi) < Q = A, y^^{ui\pk\u.i)(vi\ak.\vt) 

k=l 

(A14) 

holds for j = 1, 2, 3. Also, notice that 

J2 C^ = Pll + P22 + P 



33- 



(A15) 



Identifying C = B^, with 



R = max mm f{P,i, Pki),mm f {Pi,, Pik) } , (A16) 



in Proposition ^, we thus arrive at the bound A^ < 
(ilfljii ■ Putting all together, we have that 

^ - (l-i?)2 ~ ^ ' ^^^^> 



where in the last step we made use of Eq. (A15) 
On the other hand. 



^ A, > ^ X^tr{\u){u\pk) ■ tY{\v){v\m) = Pkh (A18) 

i i 

for all k^l =^ 1, 2, 3, and so the inequality 
4(Pii -f P22 + P33) 



Pk,i > 0, 



(A19) 



with R given by Eq. (AI6), must hold for all k, I. 

It is immediate to see that the example P in Section 
III violates it by an amount of —1/4. 
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